r 



What Is Claimed Is: 



1 LA method for using a computer system to solve a global inequality 

2 constrained optimization problem specified by a function / and a set of inequality 

3 constraints pi{x) < 0 {i=l, wherein / and pi are scalar functions of a vector 

4 X = (x;, X2, xj, ... Xrt), the method comprising: 

5 receiving a representation of the function / and the set of inequality 

6 constraints at the computer system ; 

7 storing the representation in a memory within the computer system; 

8 performing an interval inequality constrained global optimization process 

9 to compute guaranteed bounds on a globally minimum value of the function^x) 

1 0 subj ect to the set of inequality constraints; 

1 1 wherein performing the interval global optimization process involves, 

12 applying term consistency to the set of inequality 

13 constraints over a subbox X, and 

14 excluding any portion of the subbox X that is proved to be 

15 in violation of at least one member of the set of inequality 

16 constraints. 

1 2. The method of claim 1 , further comprising: 

2 linearizing the set of inequality constraints to produce a set of linear 

3 inequality constraints with interval coefficients that enclose the nonlinear 

4 constraints; 

5 preconditioning the set of linear inequality constraints through additive 

6 linear combinations to produce a preconditioned set of linear inequality 

7 constraints; 
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applying term consistency to the set of preconditioned linear inequality 
constraints over the subbox X, and 

excluding any portion of the subbox X that violates any member of the set 

of preconditioned linear inequality constraints, 

3. The method of claim 2, further comprising: 

keeping track of a least upper hound fjyar of the function^x) at a feasible 
point X wherein pi(x) < 0 (i^l, and 

including/x) < fjbar in the set of inequality constraints prior to 
linearizing the set of inequality constraints. 

4. The method of claim 2, further comprising removing from 
consideration any inequality constraints that are not violated by more than a 
specified amount for purposes of 2ipplying term consistency prior to linearizing 
the set of inequality constraints. 

5. The method of claim 1 , wherein performing the interval global 
optimization process involves: 

keeping track of a least upper bound/^Z?ar of the functionXx) at a feasible 
point x; 

removing from consideration any subbox for whichy(x) >f_bar; 
applying term consistency to the fJbar inequality /x) <f_bar over the 
subbox X; and 

excluding any portion of the subbox X that violates the f_bar inequality. 
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1 6, The method of claim 1 , wherein if the subbox X is strictly feasible 

2 (Pi(X) < 0 for all i=l, ...,n\ performing the interval global optimization process 

3 involves: 

4 determining a gradient g(x) of the function^x), wherein g(x) includes 

5 components gi(x) (i'^l, n)\ 

6 removing from consideration any subbox for which g(x) is bounded away 

7 from zero, thereby indicating that the subbox does not include an extremum of 

8 /x); and 

9 applying term consistency to each component gi(x)=0 (/= 1, of g(x)=0 

1 0 over the subbox X; and 

1 1 excluding any portion of the subbox X that violates any component of 

12 g(x)=0. 

1 7. The method of claim 1 , wherein if the subbox X is strictly feasible 

2 (p/(X) < 0 for all i=l, performing the interval global optimization process 

3 involves: 

4 determining diagonal elements (/=1, of the Hessian of the 

5 functionXx); 

6 removing from consideration any subbox for which i/„(x) a diagonal 

7 element of the Hessian over the subbox X is always negative, indicating that the 

8 function /is not convex over the subbox X and consequently does not contain a 

9 global minimum within the subbox X; 

10 applying term consistency to each inequality Hu{\) > 0 (/=1, over the 

1 1 subbox X; and 

12 excluding any portion of the subbox X that violates a Hessian inequality. 
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8. The method of claim 1, wherein if the subbox X is strictly feasible 
(pt(X) < 0 for all i=l „.,n\ performing the interval global optimization process 
involves: 

performing the Newton method, wherein performing the Newton method 
involves, 

computing the Jacobian J(x,X) of the gradient of the 
function/evaluated with respect to a point x over the subbox X, 
computing an approximate inverse B of the center of 

J(x,X), 

using the approximate inverse B to analytically determine 
the system Bg(x), wherein g(x) is the gradient of the function/x), 
and wherein g(x) includes components g,(x) (/=1, ...,«); 
applying term consistency to each component (Bg(x))^ = 0 (/-i, ,..,n) for 
each variable (/=7, ..,,n) over the subbox X; and 

excluding any portion of the subbox X that violates a component. 

9. The method of claim 1 , wherein applying term consistency 
involves: 

symbolically manipulating an equation within the computer system to 
solve for a term, g(x )\ thereby producing a modified equation g(x )) = /z(x), 
wherein the term g(x )) can be analy1:ically inverted to produce an inverse function 

substituting the subbox X into the modified equation to produce the 

equation g(X )) = h(X); 

solving for X ; = g~\h(X)); and 

intersecting X ) with the j-th element of the subbox X to produce a new 
subbox X"^; 
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12 
13 
14 



wherein the new subbox contains all solutions of the equation within 
the subbox X, and wherein the size of the new subbox X"^ is less than or equal to 
the size of the subbox X. 



1 1 0. The method of claim 1 , further comprising performing the Newton 

2 method on the John conditions. 

1 1 1 . A computer-readable storage medium storing instructions that 

2 when executed by a computer cause the computer to perform a method for using a 

3 computer system to solve a global inequality constrained optimization problem 

4 specified by a function /and a set of inequality constraints pj{x) < 0 (i=l 

5 wherein /is a scalar function of a vector x = (xj, X2, xs, ... x„), the method 

6 comprising: 



Iff ^ receiving a representation of the function/and the set of inequality 

p 8 constraints at the computer system; 

f 9 storing the representation in a memory within the computer system; 

f y 1 0 performing an interval inequality constrained global optimization process 

; J 1 1 to compute guaranteed bounds on a globally minimum value of the functionXx) 

j=f 12 subject to the set of inequaUty constraints; 

1 3 wherein performing the interval global optimization process involves, 

1 4 applying term consistency to the set of inequality 

1 5 constraints over a subbox X, and 

1 6 excluding any portion of the subbox X that is proved to be 

1 7 in violation of at least one member of the set of inequality 

1 8 constraints. 
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1 12. The computer-readable storage medium of claim 1 1 , wherein the 

2 method further comprises: 

3 linearizing the set of inequality constraints to produce a set of linear 

4 inequality constraints with interval coefficients that enclose the nonlinear 

5 constraints; 

6 preconditioning the set of linear inequality constraints through additive 

7 linear combinations to produce a preconditioned set of linear inequality 

8 constraints; 

9 applying term consistency to the set of preconditioned linear inequality 

1 0 constraints over the subbox X, and 

1 1 excluding any portion of the subbox X that violates any member of the set 

1 2 of preconditioned linear inequality constraints . 

1 13. The computer-readable storage medium of claim 12, wherein the 

2 method further comprises: 

3 keeping track of a least upper hound fjbar of the functionXx) at a feasible 

4 point X wherein ^^(x) < 0 (f =7, . . . , m) ; and 

5 including /x) < f_bar in the set of inequality constraints prior to 

6 linearizing the set of inequality constraints. 

1 14. The computer-readable storage medium of claim 1 2, wherein the 

2 method further comprises removing from consideration any inequality constraints 

3 that are not violated by more than a specified amount for purposes of applying 

4 term consistency prior to linearizing the set of inequality constraints. 

1 15. The computer-readable storage medium of claim 1 1 , wherein 

2 performing the interval global optimization process involves: 
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1 keeping track of a least upper hound fjbar of the functionXx) at a feasible 

2 point x; 

3 removing from consideration any subbox for which/x) >f_bar; 

4 applying term consistency to the fJbar inequality /(x) < fJbar over the 

5 subbox X; and 

6 excluding any portion of the subbox X that violates the fjbar inequality. 

1 1 6. The computer-readable storage medium of claim 1 1 , wherein if the 

2 subbox X is strictly feasible (pi(X) < 0 for all i=l, .,.,n\ performing the interval 

3 global optimization process involves: 

4 determining a gradient g(x) of the functionXx), wherein g(x) includes 

5 components g;(x) (/=7, ...,/?); 

6 removing from consideration any subbox for which g(x) is bounded away 

7 from zero, thereby indicating that the subbox does not include an extremum of 

8 Xx); and 

9 applying term consistency to each component gi(x)=(? {i=l, of g(x)=0 

1 0 over the subbox X; and 

1 1 excluding any portion of the subbox X that violates any component of 

12 g(x)=0. 

1 17. The computer-readable storage medium of claim 1 1 , wherein if the 

2 subbox X is strictly feasible (p,(X) < 0 for all z=7, performing the interval 

3 global optimization process involves: 

4 determining diagonal elements /f„(x) (/=1, of the Hessian of the 

5 function^x); 

6 removing from consideration any subbox for which //„(x) a diagonal 

7 element of the Hessian over the subbox X is always negative, indicating that the 
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function /is not convex over the subbox X and consequently does not contain a 
global minimum within the subbox X; 

applying term consistency to each inequality > 0 (i=l,.,.,n) over the 
subbox X; and 

excluding any portion of the subbox X that violates a Hessian inequality. 

18. The computer-readable storage medium of claim 1 1 , wherein if the 
subbox X is strictly feasible (p,(X) < 0 for all ?=i, performing the interval 
global optimization process involves: 

performing the Newton method, wherein performing the Newton method 
involves, 

computing the Jacobian J(x,X) of the gradient of the 
function / evaluated with respect to a point x over the subbox X, 
computing an approximate inverse B of the center of 

J(x,X), 

using the approximate inverse B to analytically determine 
the system Bg(x), wherein g(x) is the gradient of the function/x), 
and wherein g(x) includes components g,(x) (i^l, .,,,n); 
applying term consistency to each component (Bg(x))f = 0 (z=i, for 
each variable ,.,,n) over the subbox X; and 

excluding any portion of the subbox X that violates a component. 

1 9. The computer-readable storage medium of claim 1 1 , wherein 
applying term consistency involves: 

symbolically manipulating an equation within the computer system to 
solve for a term, g(x thereby producing a modified equation g(x )) = h{x\ 
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5 wherein the term g(x )) can be analytically inverted to produce an inverse function 

6 g'(y); 

7 substituting the subbox X into the modified equation to produce the 

8 equation g(X )) = h(X); 

9 solving for X = g~\h(X)); and 

1 0 intersecting X ) with the y-th element of the subbox X to produce a new 

11 subbox X*^; 

12 wherein the new subbox X contains all solutions of the equation within 

1 3 the subbox X, and wherein the size of the new subbox X is less than or equal to 

14 the size of the subbox X. 

1 20. The computer-readable storage medium of claim 1 1 , wherein the 

2 method further comprises performing the Newton method on the John conditions. 

1 21 . An apparatus for using a computer system to solve a global 

2 inequality constrained optimization problem specified by a function / and a set of 

3 inequality constraints pi(x) < 0 {i^L wherein /is a scalar function of a 

4 vector x = (x/, X2, xs, ... x^), the apparatus comprising: 

5 a receiving mechanism that is configured to receive a representation of the 

6 function / and the set of inequality constraints at the computer system; 

7 a memory within the computer system for storing the representation; 

8 a global optimizer that is configured to perform an interval inequality 

9 constrained global optimization process to compute guaranteed bounds on a 

1 0 globally minimum value of the function^x) subject to the set of inequality 

1 1 constraints; 

12 a term consistency mechanism within the global optimizer that is 

13 configured to, 
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apply term consistency to the set of inequality constraints 
over a subbox X, and to 

exclude any portion of the subbox X that is proved to be in 

violation of at least one member of the set of inequality constraints. 

22. The apparatus of claim 21 , further comprising: 

a linearizing mechanism that is configured to linearize the set of inequality 
constraints to produce a set of lineajr inequality constraints with interval 
coefficients that enclose the nonlinear constraints; and 

a preconditioning mechanism that is configured to precondition the set of 
linear inequality constraints through additive linear combinations to produce a 
preconditioned set of linear inequality constraints; 

wherein the term consistency mechanism is configured to, 

apply term consistency to the set of preconditioned linear 
inequality constraints over the subbox X, and to 

exclude any portion of the subbox X that violates any 
member of the set of preconditioned linear inequality constraints. 

23 . The apparatus of claim 22, wherein the global optimizer is 
configured to: 

keep track of a least upper hound fj)ar of the function/x) at a feasible 
point X wherein /?,(x) < 0 ...,m); and to 

include/x) < f_bar in the set of inequality constraints prior to linearizing 
the set of inequality constraints. 

24. The apparatus of claim 22, wherein the term consistency 
mechanism is configured to remove from consideration any inequality constraints 
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that are not violated by more than a specified amount for purposes of applying 
term consistency prior to linearizing the set of inequality constraints. 

25 . The apparatus of claim 2 1 , 
wherein the global optimizer is configured to, 

keep track of a least upper bound j^iar of the function/x) 
at a feasible point x, and to 

remove from consideration any subbox for which 
f[x) >f_ban 

wherein the term consistency mechanism is configured to, 

apply term consistency to the fjbar inequality /x) <fj>ar 

over the subbox X, and to 

exclude any portion of the subbox X that violates the/^&ar 

inequality. 

26. The apparatus of claim 2 1 , wherein if the subbox X is strictly 
feasible {p^{X) < 0 for all /-I, ,..,n): 

the global optimizer is configured to, 

determine a gradient g(x) of the function/(x), wherein g(x) 
includes components gi{\) {i^l, and to 

remove from consideration any subbox for which g(x) is 
bounded away from zero, thereby indicating that the subbox does 
not include an extremum ofy(x); and 
the term consistency mechanism is configured to, 

apply term consistency to each component gi(x)=0 
(/=1, of g(x)=0 over the subbox X, and to 
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exclude any portion of the subbox X that violates any 
component of g(x)=0'. 

21, The apparatus of claim 21, wherein if the subbox X is strictly 
feasible (p,(X) < 0 for all i^l .,.,n)\ 

the global optimizer is configured to, 

determine diagonal elements Hjiiii) (/=1, of the 
Hessian of the functiony[x), and to 

remove from consideration any subbox for which //^(x) a 
diagonal element of the Hessian over the subbox X is always 
negative, indicating that the function /is not convex over the 
subbox X and consequently does not contain a global minimum 
within the subbox X; and 
the term consistency mechanism is configured to, 

apply term consistency to each inequality ///,(x) > 0 
(i=l, over the subbox X, and to 

exclude any portion of the subbox X that violates a Hessian 
inequality. 

28. The apparatus of claim 21, wherein if the subbox X is strictly 
feasible (piiX) < 0 for all i=l, ...,n): 

the global optimizer is configured to perform the Newton method, wherein 
performing the Newton method involves, 

computing the Jacobian J(x,X) of the gradient of the 
function / evaluated with respect to a point x over the subbox X, 

computing an approximate inverse B of the center of 
J(x,X), and 
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using the approximate inverse B to analytically determine 
the system Bg(x), wherein g(x) is the gradient of the functionXx), 
and wherein g(x) includes components gj{x) {i=h and 
the term consistency mechanism is configured to, 

apply term consistency to each component (Bg(x)); = 0 
(z=i, .,.,«) for each variable over the subbox X, and to 

exclude any portion of the subbox X that violates a 
component. 

29. The apparatus of claim 21, wherein the term consistency 
mechanism is configured to: 

symbolically manipulate an equation within the computer system to solve 
for a term, g(x )\ thereby producing a modified equation g(x )) = h(x\ wherein 
the term g{x )) can be analytically inverted to produce an inverse function g ^(y); 

substitute the subbox X into the modified equation to produce the equation 

solve for X ) = g''\h{X)); and 

intersect X ) with the y-th element of the subbox X to produce a new 
subbox X"^; 

wherein the new subbox X^ contains all solutions of the equation within 
the subbox X, and wherein the size of the new subbox X is less than or equal to 
the size of the subbox X. 

30. The apparatus of claim 21 , wherein the global optimizer is 
configured to apply the Newton method to the John conditions. 
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